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Abstract 



The cross section for threshold production of the hypertriton in the re- 
action pd — > |HK + is calculated in a two-step model and compared to the 
break-up process pd — > dAK + . The latter process is shown to be dominant 
already at 2 MeV above threshold. The amplitude squared at threshold for 
the pd — > aHK + reaction is \f\ 2 = 1.0 nb/sr. 
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1 Introduction 



One way to gain information about the hyperon-nucleon interaction is to 
study bound states of the hyperon, i.e., hypernuclei. The most suitable can- 
didate for such a study would be the lightest of the hypernuclei, ^H, since the 
theoretical treatment becomes considerably easier for few-body systems (no 
collective effects that screens the perturbation caused by the hyperon). How- 
ever, the hypertriton is very weakly bound (Eb= 0.13 MeV), which makes 
the distinction between the bound and the scattering state difficult. Because 
of the weak binding, the decay properties of the hypertriton are essentially 
those of the hyperon itself. 

The purpose of the present study is to calculate the cross section for the 
process 

pd - iHK + (1) 

near threshold (Tp ab = 1126.517 MeV). In the similar reaction pd — > 3 Her?, 
the two-step model of |IJ has been successful in explaining the experimen- 
tal data. Hence, the same mechanism is employed here for the hypertriton 
production. The fundamental advantage of the two-step model is that it 
allows the exchanged momentum to be shared between the two nucleons of 
the deuteron, thereby increasing the reaction probability. In addition, the 
kinematics is miraculous in the sense that the subprocesses can be considered 
as real, which simplifies the calculations. 

The concept of a transition matrix, with a definition close to that of 
Kondratyuk et al. [0 and Komarov et al. 0, gives a somewhat different 
formulation than the one of ||. In |§, the cross section for hypertriton 
production is calculated for several angles and preferently at energies far 
from threshold, while this study is devoted to the behavior near to threshold 
and the relation to the corresponding break-up reaction. 

The index r will in this paper stand for the (A-)hypertriton. 

2 Bound final state 

For the case of a bound state between the deuteron and the lambda hyperon, 
a two step model like that of |l| is used. The Feynman diagram for this 
reaction is shown in Fig. |I[ The trick used in calculating the cross section 
for this diagram is to consider the processes of the A and C vertices to be 
almost real. This is accomplished by putting the intermediate deuteron and 
neutron on their mass shells. 
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Figure 1: Feynman diagram of the two-step model for the pd — > |HK + 
reaction, defining the various momenta. The partition of the hypertriton 
momenta is weighted by 7 = m\/(m^ + m\) = 0.373. The diagram for the 
7T° case follows upon exchanging p' and n. 



2.1 Cross section for bound final state 

the cross section in the 



With the normalization of Bjorken and Drell 
overall cm. system can be written 
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where nij and pj are the masses and momenta of the indicated particles and 
s t k is the cm. energy squared. An isospin correction factor Nj = 9/4 is 
included because of the contribution from 7r°-exchange. 

The matrix element is the integral over the Fermi momenta 



M 



d 3 k d 3 q 
(2tt) 3 (2vr) : 



(3) 



after the go an d &o integrations have been performed (giving a factor m/E w 1 
for each of the on-shell propagators when the Lorentz boost is neglected). 

The expression for M in the low energy limit and S-wave approximation 
becomes 



M = M D M C M B M A 



pi - ml + ie ' 



(4) 



3 



where the respective vertex matrices are given by 
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The parametrizations for the B and D vertices stem from vertex functions 
in the pole approximation, while for A, A and B are invariant functions of 
the total energy only and for C, Q is an invariant function of energy and 
angle. Only S-wave contributions are considered for the deuteron and the 
hypertriton. 

Summation over internal spins and polarizations gives 
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2m d 
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(9) 



where the transition matrix Gb is defined as 

■d\d 3 k #t(k)$ d (q) 



(27T)3 (p0)2 - p2 _ m 2 + le 



(10) 



The zeroth component of the virtual pion four-momentum is defined by the 
on-shell energies: p° = E d {p d )+E p {-p d )-E d {-{l^)p K -'k)-E n {\p d +'k). 
Since B is much smaller than A @, it is possible to factorize, approximately, 
the total cross section in the forward direction as 
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The cross section for the vertex processes are given by 
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3 Break-up 



Because of the small binding energy of the hypertriton, there is a large prob- 
ability that the final deuteron and hyperon are scattered instead of forming 
a bound state. This demands a calculation of the scattering cross section to 
be compared with the bound state cross section. The Feynman diagram for 
the break-up reaction is shown in Fig. [| where the kinematics is the same 
as before. 




Figure 2: Feynman diagram of the break-up reaction, pd — > dAK , defining 
the various momenta. The partition of the d and A momenta is weighted 
by 7 = mA/(md + hia) = 0.373. The diagram for the n° case follows upon 
exchanging p' and n. 



3.1 Cross section 

The kinematic part of the cross section is calculated in the low energy (non- 
relativistic) limit which gives 



da m T m p 



■ijNn fd 3 P IpkI^E \M\ 2 , (14) 



spins 



d£l 8(27r) 5 m d s T K |pd| 
where P is the relative momentum of the lambda-deuteron system and 



/ P 2 
IpkI = \/2/i K -Ad(Q - ^ — ) ( 15 ) 
V ^/Md 



is the momentum of the kaon in the final state. 

The matrix element is now given, in the low energy limit, by 



M 



(2tt) 



(16) 



where the go _ i n t e g ra tion over the neutron propagator already is performed 
and 



M = M C M B M A 



(17) 



m£ + it 

where the vertex matrices are given by Eqs (0)-©. Summing over internal 
spins and expanding the products of the Pauli matrices now gives 

-£ v /m^G S 77j v [-£p-e d cr-4, + Q-e\,]% 
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where Gs is the transition matrix element 
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(2tt) 3 / 2 (pQ) 2 - p 2 - m 2 + ie 

and the vector Q is defined by 

Q = p[^(<r-e d )] + (pXe d )[iB) + e A [B(cr-p)}. (20) 

The spinor part of Eq. ([TJ) can be separated into two orthogonal parts with 
the total spin of the deuteron-lambda system equal to 1/2 and 3/2 respec- 
tively. This is accomplished by using the projection operators 



Pi 



Pa 

2 



(21) 
(22) 
(23) 



where the two spin-1/2 projections are orthogonal (after squaring and sum- 
ming over external spins) to the spin-3/2 projection, but not to each other. 
If the B terms are neglected, the final form for the matrix element squared 
is 



6m d \Ag\ 2 (l\G s i\ 2 + l\G 
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spins 



The total cross section in the forward direction becomes 
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where the index s denotes the different spin states 1/2 and 3/2 and a s the 
corresponding fractions 1/3 and 2/3. 



4 Approximation of transition matrices 

Introducing a new parameter k 2 = (k ) 2 — m^, where k = £>°(k = q = 0), 
the pion propagator can be rewritten with the aid of the integral 

i(K+ie)r _a 

d 3 xe ip«* = ' . (26) 

r k — p^ + ie 

Substituting this for the propagator in the transition matrices Gb and Gs, s 
they become (in configuration space) 

G B = ^ />xV4(x)— e*> d (x) (27) 
47r J r 



G s , s = — Al e l P> d (x), (28) 

An J r 

where p is the pion momentum when Fermi momenta are neglected. The 
similarity of the expressions (0) and (pS| ) suggests the way to include a scat- 
tering wave function (not equal to unity) for the unbound lambda-deuteron 
system. In order to do this properly, a bound state wave function related 
to the scattering state wave function is needed. A good and fairly simple 
candidate for the bound state is the wave function of Congleton 0, which, 
unfortunately, is not easily transformed to a scattering wave function. A 
simpler alternative would be a two-term Hulthen wave function, but its be- 
haviour close to the origin is like a — br, while the Congleton function goes like 
a — cr 2 . This discrepancy could be remedied by using a modified Hulthen 
wave function (the result of using a squared Yamaguchi form factor in a 
separable potential) of the type 



ijjJr) = N- 
r 



_ e -f3 T r _ Pt + fe2 r -fi T r 

2(3 T 



(29) 



and then adjust the parameter (3 T to give the same rms radius as Congleton. 
This leads to (3 T = 23.77a r (a T = ^2fi dA E B = 0.068 fm -1 ). For the deuteron 
the S-wave part of the parametrized Paris wave function M is employed. 



4.1 Scattering wave functions 

In scattering theory it is possible to relate the scattering wave function to 
the bound state wave function. The expression for this is deduced by Faldt 
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and Wilkin |9j and is given by 



ip(r) = — 



ta{a2 + k2 k { -\k,r) 
2tt 



(30) 



k=ia 



The scattering state wave function for the spin-1/2 Ad-state has a form 
reminding of Eq. (|29|): 



^'(k,r) = — 
T 2 kr 

where 
tan 5 



(3 2 + k 2 

cos S sin kr + sin 5(cos kr — e~^ rT 1 —— — re 



-16(3 5 k(a + PY 



(31) 



(5/? 6 - 15/? 4 fc 2 - 5[3 2 k* - k 6 )(a + (3Y - ((3 2 + F) 4 ((2/? + a) 2 + (3 2 ) ' 

(32) 

For the spin-3/2 Ad-state, the same formula? are used but with another 
value for a T . The new value is estimated from the expression for the potential 
strength 

1 (3\a + (3) A 
4lrN 2 (2f3 + a) 2 + f3 2, { ' 

where N is a suitable normalization constant (independent of a). It is stated 



10|1 that the spin-3/2 case is unbound up to 1.32 times the strength pa- 



in 

rameter. Keeping f3 fixed, this gives an estimate for the strength parameter 
for a bound state (a = 0) which can be divided by 1.32 and used to calculate 
the new a' T for the unbound state. The calculations indicated above result 
in the values: 

a T = 0.068 fm" 1 
fa = 1.616 far 1 
a' T = -0.1334 far 1 . 

The wave functions for the bound state and the scattering states are plotted 
in Fig. 131. 



4.2 Final expression for the transition matrices 



The angular integrations in Eqs (]27|) and (^8|) can be performed in the S-wave 
approximation, giving (to first order), 
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5 Parametrization of vertices 



The cross sections for the vertices A and C, which turn up on the right hand 
sides of Eqs ([□]) and (|2~5"D , are parametrized empirically. 



5.1 Empirical fit of pp — > 7r + d data 



Faldt and Wilkin |1 Ij use cross section data for the reaction 7r + d — > pp and 
fit them to an exponential form. By the principle of detailed balance the 
pp — > 7r + d and 7r + d — > pp reactions have a simple relationship: 
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|Ptt| dfl 
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0.000065(^- 500)} 



(36) 



where the exponential fit is inserted in the last step. In this expression T v is 
the laboratory kinetic energy of the pion (in MeV) in the reaction 7r + d — > pp, 
giving the cross section in mb/sr. The threshold for the pd — > aHK + reaction 
corresponds to T w = 419.5 MeV, |p^| = 428 MeV/c and |p p | = 727 MeV/c, 
yielding a cross section of 53.3 /zb/sr according to Eq. (|36D . The Fermi 
momenta are neglected in this calculation. 



5.2 Empirical fit of 7r + n — ► AK + data 



There is a recent fit of the cross section for the reaction 7r + n — > AK + to a 



simple formula [12], but this does not have the correct energy dependence 
near threshold, which is crucial for the present application. Instead the 



threshold fit of Jones et al. [13 is used. This gives 



^(vr+n - AK+) = ! 
dQ K ! An 



Ap K + 



Bpl 



1 + (RpkY 



p l * b < 970 



(37) 



where A = 122 • 10" 6 fm 2 /MeV, B = 12 • 10" 9 fm 2 /MeV and R = 7.165 • 
10~ 3 (MeV) -1 (the Compton wavelength of the pion). At threshold the 
amplitude squared to use in Eqs ( |TTp and Q2"5| ) is 50.3 /zb/sr. 
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6 Results and conclusions 



The differential cross sections for the two processes pd — > iHK + and pd — > 
dAK + have been calculated in the forward direction for energies up to Tp ab = 
1137 MeV, i.e., 10 MeV above threshold. The wave functions employed 
are discussed in Sec. §. The result of the calculation is plotted in Fig. |j. 
In addition the threshold amplitude squared (|/| 2 = pi n /pft n da/dQ) for the 
reaction leading to a bound final state has been computed and found to be 

|/(pd -> A HK + )| 2 = 1.03 nb/sr. 

The cross section for the pd — > ^HK + reaction is small and in addition 
almost completely drowned by the break-up reaction pd — > dAK + . Thus, if 
one has the ambition to study a bound final state, it is of utmost importance 
to identify the hypertriton. The ratio of the cross sections for the bound and 
unbound final states in the Ad spin- 1/2 channel is given by [H 

3 



qi(dAK) ^ i / QV ^ , / 1 7W7 ^ \- 2 
a(lHK) ~ 4 U B 



.,„ -t(^) (I + a/I + Q/Eb) , (38) 



where the Q-value refers to the reaction pd — > dAK + . It is evident from this 
formula that the break-up reaction will dominate when Q is greater than Eb 
and that this will happen at small Q because of the small binding energy 
(Eb = 0.13 MeV). The ratio of the cross sections calculated from Eqs Q25D 
and (|i~l~D for the Ad spin-1/2 case agrees with Eq. (|38|) . The cross section 
for the unbound Ad spin-3/2 channel is larger than the cross section for the 
unbound Ad spin-1/2 channel, by a factor 0.7 near threshold and increasing 
to 1.9 at Tl ab = 1137 MeV. 

Some future improvements should include more accurate wave functions 



(like the one of Miyagawa et al. |i~0fl), higher energies and angular dependence. 
The deuteron wave function should be corrected for the Lorentz boost. It 
would also be valuable to have a better experimental determination of the 
7r + n — > AK + vertex. 
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Figure 3: Comparison of the Congleton hypertriton wave function (solid 
line) with the spin-1/2 (dashed line) and spin-3/2 (dotted line) Ad scattering 
wave functions. The scattering wave functions are normalized according to 
Eq. fl30D , which is approximatively valid even for k ^ ia (actually k = in 
this plot). 
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Figure 4: The differential cross sections for the reactions pd — > aHK + (solid 
line) and pd — > dAK + (dashed line) in the forward direction. 
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